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First detection of triply-deuterated methanol 
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Abstract. We report the first detection of triply-deuterated methanol, with 12 observed transitions, towards the 
low-mass protostar IRAS 16293—2422, as well as multifrequency observations of 13 CHsOH, used to derive the 
column density of the main isotopomer CH3OH. The derived fractionation ratio [CD30H]/[CH30H] averaged on a 
10" beam is 1.4%. Together with previous CH2DOH and CHD2OH observations, the present CD3OH observations 
are consistent with a formation of methanol on grain surfaces, if the atomic D/H ratio is 0.1 to 0.3 in the accreting 
gas. Such a high atomic ratio can be reached in the frame of gas-phase chemical models including all deuterated 
isotopomers of H3" . 
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1. Introduction 

Despite the relatively low elemental abundance of deu- 
terium in space (a factor of ~ 1.5xl0~ 5 less abundant 
than H; |Linsky| 1998), extremely large amounts of 
doubly-deuterated formaldehyde (D 2 CO/H 2 CO - 10%) 
have been observed in the solar-type protostar 
IRAS 16293-2422 (h ereafte r IRAS 16293,ICeccarelli et a"T1 
1998, ILoinard et aD 2000, ICeccarelli etaD 2001), initi- 
ating the search for other multiply deuterated mole- 
cules. Subsequently, doubly deuterated formaldehyde, 
doubly deuterated hydrogen sulfide and multiply 
deuterated ammonia have been observed in other pro- 
tostars and dark clouds from where protostars form 
URoueff et aO 2000. ILoinard et al.l 2001 , ICeccarelli eFaTI 
2002. lvan der Tak et al.l 2002. IUsld T 1dl 2002. IVastel et al.l 
2003). These studies have been interpreted in terms of 
two different routes for formaldehyde, hydrogen sulfide 
and ammonia deuteration: active grain chemistry fol- 
lowed by at least partial desorption into the gas for 
formaldehyde and hydrogen sulfide on the one hand and 
gas-phase chemistry for ammonia on the other hand. 
However, ammonia may also be a grain surface product, 
provided a large D/H atomic ratio in the accreting 
gas. Recently, doubly-deuterated methanol was detected 
towards IRAS 16293 l|Parise et al.l 2002). This obser- 
vation provided new constraints for chemical models. 
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The observations of the deuterated methanols CH2DOH 
and CHD2OH were both consistent with the formation 
of methanol from successive hydrogenations of CO by 
reaction with atomic H on grain surfaces, but required 
an atomic D/H ratio of 0.2 to 0.3 in the accreting gas. 
At the time of the observation of doubly-deuterated 
methanol, no gas-phase model was able to predict such a 
high atomic D/H ratio. Meantime, observations of doubly 
deuterated formaldehyde in a sample of pre-stellar cores 
showed that the degree of deuteration increases with 
increasing CO depletion IjBacmann et al.l 2002, 2003). 
This deuteration of formaldehyde in pre-stellar cores 
may occur partially in the CO-depleted gas-phase and 
partially on the surface of dust grains, followed by some 
inefficient desorption mechanism. A further spectacular 
confirmation of enhanced deuteration in CO-depleted 
gas came from the detection of abundant H 2 D + , likely 
the most abundant ion, in the prestellar core L1544 
HCaselli et al.l 2003). |Phillips and Va"stel| (2003) suggested 
that in CO-depleted gas, even the multiple deuterated 
forms of H^" may be abundant and play a role in the 
molecular deuteration enhancement. The suggestion has 
been fully confirmed by the modelling of IRoberts et al.l 
(2003), which shows that including HD^ and D^" in the 
chemical network increases dramatically the molecular 
deuteration, and allows the production of the large 
atomic D/H ratio predicted by the methanol observations 
(jParise et al.l 2002). 
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In this paper, we report the first detection of triply- 
deuterated methanol CD 3 OH in space, performed towards 
the solar-type protostar IRAS 16293. We also present a 
multifrequency observation of 13 CH30H, used to derive 
the column density of the main isotopomer CH3OH. These 
observations provide yet another stringent test to confirm 
the validity of grain surface models. 

2. Observations and results 

Using the IRAM 30-meter telescope (Pico Veleta, 
Spain), we detected the 12 CD3OH lines reported in 
Table ^ The telescope was pointed at the coordinates 
a(2000) = 16 /l 32 m 22.6 s and 5(2000) = -24°28'33.0". The 
observations were performed in April 2003. Two receivers 
were used simultaneously at 2 mm, to observe two bands 
around 156 and 160 GHz, with typical system tempera- 
tures of about 230 and 250 K respectively. These receivers 
were connected to the VESPA autocorrelator divided in 
six units. The telescope beam width is approximately 15" 
at 160 GHz. All observations were performed using the 
wobbler switching mode with an OFF position 4' from 
the source. The pointing accuracy was monitored regu- 
larly on strong continuum sources, and was found to be 
better than 3". All spectra were obtained with an integra- 
tion time of 750 minutes. The rms noise is equal to 8 mK 
(T m b) for a spectral resolution of 0.3 km s _1 . 

Observed spectra are shown in Fig. ^ The measured 
intensities, linewidths and main-beam temperatures are 
reported in Tabled The frequencies of all detected lines 
have previously been measured in the laboratory with an 
accuracy of 25 kHz (Wa lsh et al.l 1998b while the tran- 
sition strengths and energy levels were estimated from 
the published spectroscopic constants IjPredoi-Cross et al~1 
1998) using the methanol program at Ohio State. 

Some of the CD3OH lines (indicated by a star in Table 
^) are close to CH2DCN lines. In that case, the intensity 
was derived by using a two-component Gaussian fit, so 
the quoted fluxes have a further uncertainty due to the 
relative line contribution. 

Two 13 CH30H lines at 156 GHz were observed simul- 
taneously to the CD3OH lines. In addition, we analysed 
330 GHz 13 CH 3 OH observations obtained using the JCMT 
in January 2000, with an rms noise of 37 mK. The beam 
size of the JCMT is 15" at the considered frequencies, i.e. 
equivalent to the beam size of the 30-meter at 160 GHz. 
Detailed information concerning the 13 CH3 0H spectra is 
presented in Tabled 

3. Derivation of the column densities 

We derived the abundance of CD3OH using the method 
of rotational diagrams. The A and E species are consid- 
ered to be linked by ion-molecule reactions that transfer 
molecules from one species to the other. We then com- 
puted one single rotational diagram for the two species, 
presented in Fig. F5] (a). We averaged the level column 
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Fig. 2. Rotational diagram for : a) CD3OH (solid line 
T rot =85 K, dashed line: T rot = 47 K, see text), b) 
13 CH 3 OH (solid line : T rot = 208 K, dashed line : T rot = 80 
K, see text). Column densities are averaged on a 10" beam. 

densities on a 10" beam, as in IParise et all (2002), fol- 
lowing the suggestion bv lvan Dishoeck et al.l ( 1995) of en- 
hanced methanol emission in the central 10" region of 
IRAS 16293. A more recent study of the spatial distri- 
bution of CH3OH was performed bv lSchoier et al.l (2002). 
and showed evidence for an abundance jump of methanol 
of two orders of magnitude in the inner part of the enve- 
lope (< 150 AU). However, in the following we consider 
averaged abundances on a 10" beam for consistency with 
the IParise et~aTl r2002) study. 

The ground E state is estimated to lie about 4.6 K 
above the ground A state. An A state has a relative spin- 
torsional weight of 11 whereas the relative spin-torsional 
weight of an E state is 16. The partition function was 
computed from the asymmetric-top approximation : 

Z(T) = 11 X Z A (T) + 16 X exp(-4.6/T) x Z E (T) 

where Z A (T) = Z E (T) = y^fg^, with A = 3.3957 K, 
B = 0.9529 K and C = 0.9247 K, as determined by 
I Walsh et al.l H998b 

By fitting a straight line to the data in the rotational 
diagram, we derive a rotational temperature of 85 ± 28 K, 
consistent with the rotational temperature of CH 3 OH 
Ijvan Dishoeck et al.l 1995). The CD3OH column density is 
(1.4±0.9)xl0 14 cm" 2 . We also derived the CD 3 OH col- 
umn density by fixing the rotational temperature to the 
one inferred from the CH2DOH and CHD2OH molecules 
(T ro t =47± 7 K. IParise etaD 2002b The CD3OH column 
density is then (7.8 ± 2.3) x 10 13 cm" 2 . Table FJ lists the 
column densities for all deuterated methanols observed in 
IRAS 16293. 
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Fig. 1. CD3OH detected lines. The intensities are reported in main-beam brightness temperature. The label "both" 
indicates the blending of one CD3OH and one CH2DCN lines. This latter CD3OH line has not been considered in the 
population diagram analysis. 



The column density of 13 CH30H was derived using 
the same method, with the molecular parameters taken 
from the Cologne Database for Molecular Spectroscopy 
HMuller et al.l 20011. The rotational dia gram is presented 
in Fig |21 (b). The inferred rotational temperature is 
208 ±70 K and the column density averaged over a 10" 
beam is (2.6 ± 1.8) xlO 14 cm -2 . We also computed the col- 
umn density for fixed rotational temperatures of 50 K 
and 80 K. The inferred value is (1.4 ± 0.6) x 10 14 cm~ 2 , 
independent of the temperature in this range. Using the 
12 C/ 13 C ratio of 70 derived by |Boogert etaTj (2000), we 
derive a column density of (9.8 ± 4.2) x 10 15 cm~ 2 for 
CH3OH. This value is nearly 3 times higher than the col- 
umn density used bv IParise et al.l (2002), inferred from 
CH3OH observations from Ivan Dishoeck et al.l (1985). 



The fractionation ratios, relative to this new estimate of 
the CH3OH column density, arc reported in Table [5] for 
all deuterated isotopomers of methanol. 



4. Discussion and conclusions 

The main result of this Letter is the first detection of 
triply-deuterated methanol in space, with 12 detected 
transitions. This discovery follows the detection of doubly- 
deuterated as well as singly-deuterated isotopomers to- 
wards the same object l|Parise et al.l 20021. Observations 
of multiple isotopomers of methanol represent a powerful 
constraining tool for chemical processes that lead to such 
a high deuteration. 
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Table 1. Main-beam intensities 1 , peak temperatures 2 and widths for the observed CD 3 OH and 13 CH 3 OH transitions 3 . 



Frequency 


Transition 






J T mb dv 


T m b 


Av 


GHz 




Debye 2 


K 


K km s^ 1 


mK 


km s _1 


CD3OH 




4:1-01 -Cj2 


2.94 


21.5 


1 09 4- fl t)9^ 


1 


1.8 ± 0.4 


156.239295 


4 2 -3 2 A- 


2.44 


42.0 


0.032 ± 0.012 


26 


1.2 ± 0.5 


156.242613 


4 -3o A+ 


3.13 


18.8 


0.059 ± 0.018 


31 


1.8 ± 0.5 


156.253079 


4 2 -3 2 A+ 


2.44 


42.0 


0.056 ± 0.015 


38 


1.4 ± 0.4 


156.260737* 


43-33 E 2 


1.39 


55.5 


0.061 ± 0.018 


46 


1.3 ± 0.4 


156.262936 


4 3 -3 3 Ei 


1.37 


46.9 


0.044 ± 0.011 


43 


1.0 ± 0.2 


156.275238 


4i-3i Ei 


2.93 


33.1 


0.034 ± 0.021 


26 


1.2 ± 1.2 


156.285288* 


4 2 -3 2 E 2 


2.35 


36.3 


0.065 ± 0.021 


54 


1.1 ± 0.3 


156.581519 


8i-7o E 2 


4.15 


70.2 


0.050 ± 0.010 


46 


1.0 ± 0.2 


160.640122 


2o-2i E 2 


2.38 


15.95 


0.090 ± 0.019 


26 


3.2 ± 0.5 


160.718291 


6 2 -5i Ei 


2.22 


50.2 


0.039 ± 0.009 


29 


1.3 ± 0.2 


160.753934 


lo-li E 2 


1.45 


12.2 


0.038 ± 0.009 


35 


1.0 ± 0.2 


1J CH 3 OH 


156.299374 


5o5-5_15 


0.697 


47.1 


0.19±0.08 


65 


3.0 ± 0.5 


160.507694 


2l2-3o 3 


0.300 


21.3 


0.08±0.02 


29 


2.4 ± 0.6 


330.194042 


7-17-6-16 


5.55 


69.0 


0.51±0.16 


120 


3.8 ± 0.6 


330.252798 


7o7"6o6 


5.66 


63.4 


0.43±0.15 


100 


4.0 ± 1.2 


330.265233 


7-61-6-60 


1.50 


253.4 


0.15±0.10 


80 


1.7 ± 1.0 


330.277270 


761-660 


1.50 


258.1 


0.12±0.08 


60 


2.0 ± 0.7 


330.277270 


76 2 -6ei 


1.50 


258.1 


0.12±0.08 


60 


2.0 ± 0.7 


330.319110 


75 2 "651 


2.78 


202.0 


0.45±0.21 


90 


4.8 ± 0.9 


330.319110 


753"65 2 


2.78 


202.0 


0.45±0.21 


90 


4.8 ± 0.9 


330.342534 


744-643 


3.81 


144.2 


0.07±0.05 


50 


1.3 ± 6.7 


330.342534 


743"64 2 


3.81 


144.2 


0.07±0.05 


50 


1.3 ± 6.7 


330.408395 


734-633 


4.61 


111.4 


0.91±0.19 


200 


4.2 ± 0.7 


330.442421 


7l6"6l5 


5.69 


84.49 


0.42±0.12 


140 


2.8 ± 0.9 


330.535822 


7 2 5"6 2 4 


5.14 


85.80 


0.25±0.11 


40 


5.6 ± 1.2 


330.535890 


7- 2 6"6_ 2 5 


5.20 


89.45 


0.25±0.11 


40 


5.6 ± 1.2 



x The fluxes were derived using Gaussian fits, and the uncertainty given is y/ tj 2 tai + o\ al where a s tat is the statistical error and 
a ca i the calibration uncertainty (15%). 2 The noise rms is 8 mK for the CD3OH data and 37 mK for the 13 CH 3 OH data. 3 A 
star following the frequency indicates that the line is close to a CH 2 DCN line and was fitted by a two-component Gaussian fit 
(see text). 



Table 2. Derived column densities and fractionation 
ratios relative to CH3OH for deuterated methanols in 
IRAS 16293 



Molecule 
CD3OH 

CHD2OH 6 
CH 2 DOH b 
CH 3 OD b 



T rot (K) 
85 ± 28 
47 ± T 

47 ±7 

48 ± 3 
20 ± 4 



N (cm" 2 ) 
(1.4 ± 0.9) xlO 14 
(7.8 ± 2.3) xlO 13 
(6.0 ± 2.2) xlO 14 
(3.0 ± 0.6) xlO 15 
(1.5 ± 0.7) xlO 1 4 



fractionation 
1.4 ± 1.4 % 
0.8 ± 0.6 % 

6 ± 5 % 
30 ± 20 % 
2 ± 1 % 



"fixed temperature, see text. ^Observed in Parise et al. 2002 

It is interesting to compare these observations to 
the predictions of the simple grain chemistry scheme 
of |Rodgers fc Charnley| (2002). If the D atoms are ran- 
domly distributed in the methanol isotopomers (i.e. this 
scheme does not consider any activation barrier for 
the reactions but rather assumes that all reactions are 
equiprobable) , the fractionation ratios R of each iso- 
topomer relative to CH3OH should scale as follows: 
R(CH 3 OD) = a, R(CH 2 DOH)=3a, R(CHD 2 OH) = 3a 2 



and R(CD 3 OH) = a 3 , where a is the accreting atomic 
D over H ratio. The three independent observations of 
CH 2 DOH, CHD 2 OH and CD 3 OH are consistent within 
the error bars with a value of 0.1—0.2 for the D over H 
accretion rate. Accounting for the different mass of the 
atoms, this ratio corresponds to an abundance ratio in the 
gas-phase of D/H = \/2x (0.1—0.2) = 0.15-0.3. However, 
this simple scheme fails to explain the observed low abun- 
dance of CH 3 OD. 

More accurate grain chemical models accounting for 
different activation barriers for the reactions have been 
developed in the last few years. We compare in the fol- 
lowing our observations with the model developed by 
IStantcheva fc Herbstl (2003). This model is based on the 
direct solution of the master equation and therefore gives 
essentially the same predictions as the Monte Carlo mod- 
els described by ICaselli et al.l (2002) or |Charnley et al.| 
(1997). Fig. |21 shows predictions for fractionation ra- 
tios of deuterated isotopomers of methanol relative to 
CH 3 OH versus the atomic D/H ratio in the accreting gas 
<|Stantcheva fc Herbstl 2003) when the various isotopomers 
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Fig. 3. The deuterium fractionation ratio for the various 
deuterated isotopomers of methanol is plotted against the 
abundance ratio of deuterium to hydrogen atoms in the 
gas phase. Dashed lines: model results of Stantcheva et al. 
2003. Thick lines: observations of IRAS 16293. 

are formed by active grain chemistry. In the limit of low 
temperature (10 K), this model essentially gives the ra- 
tios corresponding to a random distribution of deuterium 
atoms. Observed fractionation ratios with their error bars 
have been overlaid on each curve, allowing the derivation 
of the required atomic D/H ratio in the gas-phase. The 
CD3OH, CHD 2 OH and CH 2 DOH observations are con- 
sistent with a formation on grain surfaces with an atomic 
D/H abundance ratio of 0.1—0.2. Such a high atomic frac- 
tionation ratio in the gas phase is predicted by the recent 
gas-phase model of IRoberts et al.l f2003). which involves 
not only H2D+ but also D2H + and Djj~ as precursors for 
deuterium fractionation, when the density of gas is very 
high and heavy species such as CO are strongly depleted. 

As can be seen in Fig. [21 CH3OD appears to be under- 
deuterated when compared with the grain chemical pre- 
dictions. It is possible that the CH3OD fractionation may 
be affected in the warm gas; e.g., this isotopomer may 
be preferentially converted into CH 3 OH when released in 
the gas-phase by protonation reactions followed by dis- 
sociative recombination with an electron ( | Charnley et al.| 
1997, EaFise..eFar 2002) : 



CH3OD + Ho 



H 2 
D. 

HD+ and D+ 



CH3ODH+ + 
CH3OHD+ + e- -► CH3OH + 
The corresponding reactions with H2D + , 
come of little importance in view of their low abundance 
in the warm gas of the hot core. This hypothesis, which as- 
sumes that protonation reactions attack the oxygen end of 
the methanol only (Osamura, Roberts & Herbst, in prep), 
could be tested by observing the CH2DOD isotopomer. 
This observation may be difficult due to the expected low 
intensity of the lines. 

Acknowledgements. E. Herbst acknowledges the support of the 
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